Light pulse in A-type cold atomic gases 
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We investigate the behavior of the light pulse in A-type cold atomic gases with two counter- 
propagating control lights with equal strength by directly simulating the dynamic equations and 
exploring the dispersion relation. Our analysis shows that, depending on the length Lq of the 
stored wave packet and the decay rate 7 of ground-spin coherence, the recreated light can behave 
differently. For long Lo and/or large 7, a stationary light pulse is produced, while two propagating 
light pulses appear for short Lo and/or small 7. In the 7^0 limit, the light always splits into two 
propagating pulses for sufflciently long time. This scenario agrees with a recent experiment [Y.-W. 
Lin, et al, Phys. Rev. Lett. 102, 213601(2009)] where two propagating light pulses are generated 
in laser-cooled cold atomic ensembles. 

PACS numbers: 32.80.Qk 42.50.Gy 
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I. INTRODUCTION 

Quantum information transfer between light and 
atomic ensembles has attracted much attention recently. 
In particular, electromagnetically induced transparency 
(EIT) [1, 2], a robust technique that renders a reso- 
nant opaque medium transparent by means of destruc- 
tive quantum interference, has been exploited to realize 
the storage and retrieval of light pulses in atomic ensem- 
bles [3-6]. In the storage process, a weak probe light 
pulse carrying quantum information and a strong cou- 
pling light are applied to an optically thick atomic en- 
semble. The probe light is then gradually converted into 
a ground-spin coherence as the coupling light is adiabat- 
ically switched off, and, as a result, the quantum infor- 
mation is stored in the atomic ensemble. The reading 
process is almost reverse to the storage: the control light 
is adiabatically turned on, and accordingly a new light 
pulse is created and propagates out of the atomic ensem- 
ble. In this way, the quantum information can be stored 
and read out without loss in principle. If two counter- 
propagating control light pulses with equal strength are 
adiabatically switched on [7], the retrieved light pulse 
will not propagate out of the atomic ensemble. Instead, 
it stops in the media and forms a stationary light pulse, 
as experimentally demonstrated by Bajcsy et al. [8]. On 
this basis, together with cold-atom techniques, many ap- 
plications are proposed, which includes simulation of the 
dynamics of massive Schrodinger particles [9] , of Dirac 
particles [10], and of the strong correlated Bose system 
confined in a hollow core fiber etc. [11]. 

The dynamic equations describing the behavior of the 
recreated light in A-type atoms contain infinite-order 
terms. In conventional theoretical treatment, secular 
approximation is used and only the zcroth-order coef- 
ficient of the ground-spin/optical coherence is kept [12]. 
In hot atomic systems like in Ref. [8], such a treatment 
is reasonable, since the higher-order terms decay very 
fast due to the random atomic motions and collisions. 



However, thermal fluctuations are strongly suppressed 
in cold atoms, and thus higher-order terms decay much 
slower and should be considered. The light generation 
in cold atomic ensembles has been theoretically studied 
beyond the secular approximation by Hansen et al. [13]. 
They came to the conclusion that, when the decay rate 
of the ground-spin/optical coherence is zero, the gener- 
ated light is a pure stationary light pulse-i.e., a station- 
ary light without photon loss. Nikoghosyan et al. [14] 
took into account the relaxation of the upper state and 
demonstrated that, under the slow-light condition, the 
generated light is a stationary light pulse but with some 
photon loss. Nevertheless, a recent experiment demon- 
strated [15] that, in laser-cooled cold atomic ensembles, 
the retrieved light pulse is not stationary but splits into 
two propagating wave packets. The authors gave a sim- 
ple model which only involves the zeroth-order and first- 
order coefficients. Since there is no obvious reason why 
the cutoff should take place at the first-order term, our 
original minor motivation was to deal with the dynamic 
equations to a higher order. It turns out that, in the 
zero-decay limit and under adiabatic approximation, the 
dynamic equations can be analytically treated to any or- 
der as one wishes (see Appendix for details). Indeed, to 
the first order, one obtains two counter-propagating light 
pulses, consistent with Ref. [15]. However, as higher- 
order terms gradually come in, the relative group ve- 
locity of the two light pulses decreases and vanishes as 
\/2£ + 1/ i, where £ is the highest order of term in the cal- 
culation. It seems to confirm the results in Refs. [13, 14]. 
Given this discrepancy and the potential important ap- 
plications of the stationary light pulse in cold atoms, a 
careful and systematic study seems desirable. 

In this paper, we directly simulate the dynamic equa- 
tions (given in Sec. II) to avoid further approximation. 
For a given set of parameters, a series of simulations is 
performed with cutoff at different order £ and the re- 
sult are extrapolated to the infinite-order limit {£ 00). 
These results are presented in Sec. HI. Section IV pro- 
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vides a qualitative understanding from the numerical cal- 
culation of the dispersion relation, which is obtained from 
the Fourier transformation of the dynamic equations. A 
brief discussion is given in Sec. V. 



II. DYNAMIC EQUATIONS FOR THE 
ATOM-LIGHT SYSTEM 



Let us consider an ensemble of A-type atoms aligning 
along a certain direction (say z), which is horizontal in 
Fig. 1. These atoms interact with with a weak probe 
light Ep and a strong control light , treated as quan- 
tum and classical light, respectively. 
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FIG. 1: (Color online) Sketch of the interaction between A- 
type atoms and control light JlJ and probe light E^. 



Under the singlc-modc approximation, the interaction 
Hamiltonian in the rotating frame reads [4] 



N f ^ 

% = -— I dzhga3i{z,t)Ep{z,t)+h(T32{z,t)V,c{z,t)+h.c., 

where g is the coupling constant, aij{z,t) := 
J dzmCfij5{z — Zm) is the continuous atomic operators, 
with a™ := (j| the spin flip operator of the mth 
atom, N is the atomic number, L is the length of the 
atomic ensemble, Ep{z,t)=Ep{z,t)e~^'^''^ is the electric 
field of the probe hght, and ^c{z, t)=flc{z, t)e~^'^''* is the 
Rabi frequency of the control field. For simplicity, we as- 
sume the two ground states are degenerate, and the probe 
light and control light are on resonance Wp « Wc = '^31- 
With slowly varying atomic operators 

ai3{z,t)^ai3{z,t)e''^''' , 

<Ti2(^,i)-2^i2(^,i)e"^-*-""=*, 

the Langevin equations governing the atomic dynamics 
read as [16] 



da- 



is 



dt 

dai2 
dt 



= -r6-i3 + igEp + iQc^i2 + Fi3 , (1) 
= iQlau + #12 , (2) 



where we have set an = 1, (T33 = a23 = 0. This approx- 
imation is appropriate since the probe field is very weak 
and all the atoms are initially prepared in F is the 
decay of the optical coherence, and Fij is the Langevin 
force. 

The counter-propagating control light can be described 
by Qciz,t) = Q+e'-''''^ + Oje"*'"''^, where the light is as- 
sumed to be homogeneous. The probe light can also be 
decomposed into two counter-propagating components as 



Ep{z,t) = E+iz,t)e'''^' + E;iz,t)e 



-ikcZ 



(3) 



Following the standard procedure [12-14], we define 
ground-spin coherence as S" = y/Na2i and optical co- 
herence as P = •\//V<5'3i, and expand them as 



,2nikcZ 



n— — 00 
n— 00 



p- E p^ 



2n+ie 



(2ri+l)ifcc2 



(4) 
(5) 



Inserting Eqs.(3)-(5) into Eqs.(l) and (2) and assuming 
= r^c, wc obtain a set of dynamic equations as 



2n+l 



dt 



dS2n 

dt 



-(r + 72n+l)P2«+l 



+ igVNEp'2n+l + i^c{S2n + S2{n+1)) (6) 



— ~l2nS2n + if^c(^2ri-l + -f2n+l) i 



(7) 



where i?p,±i = E^ , Ep^2n+i{n 7^ 0,-1) = 0, and 7„ 
represents the decay of nth-order coefficient. We have 
also neglected the Langevin force terms since they do 
not play a role in the long-time behavior of the light 
pulse. The dynamics of the probe light is governed by 
the Maxwell equations 



d_El 
dt 

dE- 
dt 



d_El 

' dz 

dEp 

^ dz 



= igVNPi , 
= ig^/NP-i 



(8) 
(9) 



For warm atomic vapors, the random motions and col- 
lisions of atoms result in a very rapid decay of spatial 
coherence. Effectively, one has 70 = and 7„ > for 
n 7^ 0. In this case, multiple components S2n{n 7^ 0) of 
ground-spin coherence are suppressed, the same applies 
to P2n+i{n 7^ —1,0). Thus, these terms can be neglected, 
and the probe light forms a stationary light pulse [12]. 

In contrast, in a deep optical lattice, where atoms are 
fixed at the lattice sites [19], the decays of the higher- 
order coefficients can be ignored. In other words, one 
has 7n = for any n. Thus, the multiple components 
can be populated and preserve their coherence [12], and 
secular approximation is no longer valid. Nevertheless, 
after adiabatic elimination, one can analytically solve the 
Eqs.(6)-(9) with 7„ = (see Appendix for details). The 
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FIG. 2: (Color online) Group velocity Vg of the recreated 
forward (backward)-propagating light pulse in the 7n = 

limit. The unit of v„ is 



The red straight line is from the 
approximate analytic calculation (24) (25) in Appendix, while 
the black dots are obtained by the direct numerical simulation 
of Eqs.(6)-(9) in Sec. III. The error margins of the data points 
are ±0.006. 



group velocity Vg of the forward (backward)-propagating 
light pulse is shown in Fig. 2. One finds that: 1), the ap- 
proximate solution with cutoff at finite £ yields a non-zero 
Vg, and 2), Vg{i) reaches the maximum value at £ = 1, 
and then vanishes with ^/2£ + l/£. Since the real sys- 
tem corresponds to the infinite-order limit, one "seems" 
to conclude that the recreated light forms a stationary 
pulse. 

For cold atomic systems where temperature is low but 
nonzero, one has 70 = and 7n 7^ for n ^ 0. In this 
case, we set 7„ = |n|ar, with a the decay constant. This 
decay model can well describe the decays of the coef- 
ficients of the ground-spin/optical coherence in various 
cold atomic systems. For instance, in the laser-cooled 
cold atomic ensembles, the higher-order coefficients have 
a phase grating of e™*'''=^ across the atomic gases, and 
thus will decay due to atomic random motion. The de- 
cay rate can be estimated by the time needed for the 
atoms moving across one wavelength of the phase grat- 
ing 7„ ~ = [17]. In the Bose condensation, 
the higher-order coefficients can be regarded as a particle 
excitation with momentum |n|fcc?i. At average, they move 
out of the atomic gases after a time of L{ l"l^'^^ )-i [ig]^ 
and the decay rate can be approximated by 7„ ~ ^""jnl^ ■ 



III. NUMERICAL SOLUTION 

To check the vafidity of the approximate solution in 
Fig. 2 and further find the dynamics of the recreated light 
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FIG. 3: (Color online) Light intensity -I- as a func- 
tion of time and position z for 7™ = with different length 
io- Figures (a)-(f) represents Lo/lats = 5,10,20,30,40, and 
50, respectively. The position is in unit of labs and the time 
is in unit of l/F. Strong light is shown in red bright color, 
while the back ground is in blue. 



pulse in cold atomic systems (7^ = |7i|ar) for which the 
approximate treatment is unavailable, we directly simu- 
late Eqs.(6)-(9). Naturally, a cutoff takes place at finite £ 
and accordingly b+M equations are involved in each sim- 
ulation. The result for real systems {£ — >■ 00) is obtained 
from the extrapolation of simulations for finite £. 

We first consider the zero-decay limit (7„ — 0). The 
initial condition is taken such that: 1), only the zeroth 
component So{z,t = 0) of the ground-spin coherence is 
nonzero while all other components S2n{z,0) are zero; 
So assumes a Gaussian shape 5*0(2;, 0) = e"^^/^"^ with 
Lq the length of the wave packet; 2), all components of 
the optical coherence are zero P2n-\-i{z ^ 0) — O5 and 3), no 
probe light exists at the beginning £"+ (z, 0) = E~ {z,0) = 
0. Further, the wave-packet length is set at Lq = blabs, 
with lahs = the absorption length. To be in the slow- 
light regime, we chose the parameters to be fie = O.GQF, 
g^N = 138r^. The equations are directly solved by Lax- 
Friedrichs method with sufficiently small step. The sim- 
ulation is up to £ = 100, and the group velocity Vg of the 
recreated forward (backward)-propagating light pulse is 
measured. The results are shown in Fig. 2. As the ap- 



proximate analysis, Vg{£) reaches its maximum at £ = 1 
and then starts to decrease. Nevertheless, the decrease 
of Vg becomes slower and slower after £ w 14 and eventu- 
ally stays unchanged at Vg = 0.47 ± 0.006. From Fig. 2, 
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FIG. 4: (Color online) Light intensity |£d|^ for Lo = 

Slabs with different decay constant a. Figures (a)-(f) are for 
a = 0.2, 0.02, 0.01, 0.005, 0.001, and 0, respectively. The posi- 
tion is in unit of lahs and the time is in unit of l/F. 



FIG. 5: (Color online) Light intensity \Esf+ \Edf for 
a — 0.001 with different length Lq. Figures (a)-(f) are for 
Lo/labs = 50,40,30,20,10, and 5, respectively. The position 
is in unit of labs and the time is in unit of l/F. 



it looks rather secure to conclude that the group velocity 
Vg takes a finite value for ^ — >■ oo, suggesting that the 
recreated light splits into two counter-propagating light 
pulses. This defies the earlier approximate solution that 
the light forms a stationary light pulse. 

We then consider the effect of the wave-packet length 
Lq for the 7„ = case and set Lo/labs = 5, 10, 20, 30, 40, 
and 50. We define the sum mode Eg = + E~ and 
difference mode Ed = E^ — E^ , and measure the inten- 
sity of the light pulses |£'sp-|- \Ed\'^ as a function of time 
and position z. The results for i = 30, where the simula- 
tion already reaches the steady state, are shown in Fig. 3. 
One observes that, for all the cases, the generated light 
always splits into two counter-propagating light pulses. 
Nevertheless, as grows, the time for the occurrence of 
splitting becomes longer and longer. In experiments that 
are finished within short times, one may not be able to 
observe such a splitting. 

Next, we study the decay model for cold atomic sys- 
tems (7„ = |n|ar). Shown in Figs. 4 (a)-(f) arc the dy- 
namics for Lo = Uabs with a = 0.2, 0.02, 0.01, 0.005, 0.001 
and 0. When the decay rate is large-i.e., a is large, the 
generated light forms a stationary light pulse with dissi- 
pating; sec Fig.4(a)-(c). In contrast, for sufficient small 
a (Fig.4(d)-(f)), two counter-propagating light pulses ap- 



pear. 

Figures 5 (a)-(f) display the dynamics of the generated 
light for a = 0.001 with Lo/labs = 50, 40, 30, 20, 10, and 
5. It can be seen that the generated light is stationary 
for long wave length Lq while splits into two light pulses 
for small Lq. Wc should mention that, however, the dy- 
namics shown in Fig. 5 is up to time lOO/F, much shorter 
than 500/r in Fig.3(c)-(f). To see the behavior of the 
light pulse for longer time, we performed the calculation 
for Lo/labs = 50 up to the time 500/r and did not ob- 
serve any splitting. 

We also calculated the total light strength / = 
/ (l-EsP + \Ed\'^)dz remaining in region {— 3io, 3Lo} as a 
function of time. The results for Lq = blabs and different 
a are shown in Fig. 6. Indeed, as a decreases, the remain- 
ing light strength / becomes weaker and weaker, reflect- 
ing that the loss of photons becomes more and more se- 
rious due to the propagating of the splitting light pulses. 



IV. DISPERSION RELATION 

The earlier simulation yields directly observable phe- 
nomena. In this section, we aim to provide a qualitative 
understanding by exploring the associated dispersion re- 
lation. 
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FIG. 6: (Color online) Remaining light strength I as a func- 
tion of time for Lq = 5labs with a = 0.2, a — 0.02, a — 0.01, 
a = 0.005, a = 0.001, and a — 0. The unit of / is arbitrary 
and the unit of time is l/F. 



Note that Eqs.(6)-(9) are linear, the Fourier transfor- 
mation of these equations leads to 



(r + 7|2 



- iUj)P2n+l = igVNEp^2n+l 

+ iflc{S2n + S2(n+1)) , (10) 
(7|2ri| - i'^)S2n = i^ciP2n-l + P2n+l) , (H) 



ickEp 



oE- - ickE:. 



: igVNP^i . 



(12) 
(13) 



After some tedious calculations, similar as those in Ap- 
pendix, we obtain the dispersion relation 



c 



\ 



r 



- -) (14) 



where we have introduced two effective decay parameters 
Psiiv) and Foiuj) 

r.iuj) = r + 71 - icj + + n , 

Fdiuj) = r + 71 - + 7?. , 
with parameter TZ recursively expressed as 

7^:= 



72 - + 



In the dispersion relation (14), the real and the imaginary 
parts of momentum fc, |Re(fc)| and |Im(fc)|, qualitatively 
characterize the effects of dispersion and of dissipation, 
respectively. 
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FIG. 7: (Color onhne) Numerical results for dispersion re- 
lation (14) for Lo = and different decay constants a. 
The vertical axis represents the real and the imaginary part 
of momentum k: K := (c/F)Im(fc) (the red thin curve) and 
(c/F)Re(fc) (the black thick curve). The insets display the 
dispersion relation in a much larger region. 



For a given frequency uj, we numerically calculate 
Re(A:) and Im(fc) by taking a cutoff at £ for TZ, and then 
extrapolate the calculations to £ — > oo. The parame- 
ters g^N = 138r^ and fie = O.GOF are the same as in 
Sec. II, and the decay model 7„ = |n|aF is used. Fig- 
ure 7 displays the results for Lq = blabs at £ = 1000, 
where TZ already reaches its steady value. In Fig. 7(a) 
where a — 0.2 is comparable to 1, the black thick (Re(fc)) 
and the red thin (Im(A:)) curves almost overlap in the 
vicinity of original point (|A:| = 0,w = 0). As w in- 
creases, they gradually separate from each other. The 
red thin curve increases faster than the black thick one- 
i.e., |Im(fc(w))|>|Re(fc(w))|. This implies that dissipation 
dominates over dispersion. When a becomes smaller and 
smaller, the two curves separate at smaller frequency uj. 
Furthermore, one has |Im(A:(a;))|<|Re(/c(a;))|. Namely, 
the effect of dispersion becomes more important than 
that of dissipation. For a < 0.01 (Fig.7(d)-7(f)), the 
two curves significantly separate from each other as long 
as |fc| deviates from zero. This means that the light pulse 
has a nonzero group velocity and can propagate out of 
the atomic ensembles. In addition, one observes that the 
black curve grows less and less rapidly as a decreases, 
reflecting that the effect of dissipation becomes weaker 
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and weaker. The insets show that, for very large |fc|, the 
dissipation always play the major role. 

In short, we argue that the phenomena in Sec. Ill 
can be qualitatively understood from the competition be- 
tween the effects of dispersion and of dissipation. For 
a « 1, the dissipation dominates over the dispersion, and 
a stationary light pulse is generated. This applies to ther- 
mal atomic gases. For a ^ 1, the dispersion wins as long 
as |fc| > 0. If the length of the stored wave packet Lq 
is sufficiently long-i.e., |fc| « 0, the recreated light forms 
a stationary pulse; otherwise, it splits into two counter- 
propagating light pulses. In the a = limit and for a 
finite length wave packet, the generated light will always 
splits into two propagating light pulses for sufficiently 
long time, since the dispersion is always dominant over 
dissipation. 



ination, as illustrated below. Assuming that the charac- 
teristic interaction time T is long compared to the upper 
level rclaxation-i.e., ^ T, wc can adiabatically elimi- 



nate in Eqs. (6) and (7) and obtain 

rP2n+l = igEp^2n+l + i^c{S2n + '5'2(n+l)) , (15) 
^-zr!,(P2„_i+P2n+l). (16) 



Differentiating both sides of Eq. (15) with respect to t 
and substituting Eq. (16) yield 

dP2n+l . dEp .2n+l 

r — - — = ig- 



dt 



dt 



-nl[P2n-l + 2P2«+1 + P2«+3) . 



(17) 



V. DISCUSSION 

In summary, using direct simulation of the dynamic 
equations for A-type atomic systems, we find that both 
the decay rate 7 and the length Lq of the stored wave 
packet play an important role in determining the behav- 
ior of the new light generated by two counter-propagating 
control lights with equal strength. The numerical simu- 
lation of the 7 = limit defies the approximate analyt- 
ical solution. This means that the adiabatic-elimination 
treatment demonstrated in Appendix is invalid. For cold 
atomic systems, our calculations suggest that the recre- 
ated light forms a stationary pulse for large Lq and/or 
7 while splits into two counter-propagating light pulses 
for small Lq and/or 7. This scenario agrees well with 
the recent experiment. A qualitative understanding is 
given from the aspect of the dispersion relation. We ex- 
pect that our systematic calculation shall provide useful 
information for future experiments. 



Once more, the derivative terms can be eliminated by 
adiabatic treatment, and the equations become 

Ps.2n+i = (-l)"((2n + l)P,,i - n^^) , (18) 

P..2n+i = (-l)"(P.,i-ng^). (19) 

where we have introduced s mode and d mode 

^!s,2n+l = P2n+1 + P^{2n + l) -i (20) 
PcL2n+l = P2n+1 — P-{2n+l) ■ (21) 

Neglecting terms Ps,2n-i-i and Pd,2n+i for n > i, one 
has 



ig dEs 



2£+inl dt ' 
, ig dEd 



Pd.i = £ 



m dt 



(22) 
(23) 



Acknowledgments 

We are very grateful to Ite A.Yu and Tao Xiong 
for helpful discussions. This work is supported by 
the NNSFC, the NNSFC of Anhui (under Grant No. 
090416224), the CAS, and the National Fundamental Re- 
search Program (under Grant No. 2006CB921900). 



Making use of the Maxwell equations and the initial 
condition Sn = for n 7^ 0, we obtain the propagating 
solution 

^^^^^^j^ ^ O.cSo{z - cot,Q) ^ Q.aSo{z + cot,Q) ^ ^^4) 
9 9 ' 

' CO g g ' 



Appendix 

For 7„ = 0, Eqs. (6) and (7) can analytically solved 
under additional approximation-i.e., the adiabatic elim- 



is 



Where the group velocity of the splitting wave packet 
Co = + + In low the group- 

'2f+isi|; 



velocity limit, cg « c- 



eg^N 
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